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(Dated: October 21,2007) 

We report results of a search for supersymmetry (SUSY) with gauge-mediated symmetry breaking 
in di-photon events collected by the DO experiment at the Fermilab Tevatron Collider in 2002-2006. 

In 1.1 fb-i of data, wc find no significant excess beyond the background expected from the standard 
model and set the most stringent lower limits to date for a standard benchmark model on the lightest 
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neutralino and chargino masses of 125 GeV and 229 GeV, respectively, at 95% confidence. 



PACS numbers: 14.80.Ly, 12.60.Jv, 13.85.Rm 

Low-scale SUSY is one of the most promising solutions 
to the hierarchy problem associated with the intrinsic 
disparity between the electroweak and Planck scales. It 
postulates that for each known particle there exists a su- 
perpartner, thereby stabilizing the radiative corrections 
to the Higgs boson mass. Bosons have fcrmion super- 
partners, and vice versa. None of the superpartners have 
yet been observed, and superpartner masses must there- 
fore be much larger than those of their partners, i.e., 
SUSY is a broken symmetry. Experimental signatures 
of supersymmetry are determined through the manner 
and scale of SUSY breaking. In models with gauge- 
mediated supersymmetry breaking (GMSB) P, IJ , it is 
achieved through the introduction of new chiral super- 
multiplets, called messengers that couple to the ultimate 
source of supersymmetry breaking and to the SUSY par- 
ticles. At colliders, assuming i?-parity conservation 3|, 
superpartners are produced in pairs {xtxi ^^^d xfx2 
production dominates in most cases) and decay to the 
standard model particles and next-to-lightest SUSY par- 
ticle (NLSP), which can be either a neutralino or a slep- 
ton. In the former case, which is considered in this 
note, the NLSP decays into a photon and a gravitino 
(the lightest superpartner in GMSB SUSY models, with 
mass less than « 1 keV). The gravitino is stable, and es- 
capes detection, creating an apparent imbalance in trans- 
verse momentum {^t) in the event. GMSB SUSY final 
states are therefore characterized by two energetic pho- 
tons and large missing transverse momentum. The dif- 
ferences in event kinematics between particular GMSB 
SUSY models result in slightly different experimental 
sensitivities [1], and to obtain a quantitative measure 
of limits on SUSY we consider a model referred to as 
"Snowmass Slope SPS 8" This model has only a 

single dimensioned parameter: an energy scale A that 
determines the effective scale of SUSY breaking. The 
minimal GMSB parameters correspond to a messenger 
mass Mm = 2 A, the number of messengers A''5 = 1, the 
ratio of the vacuum expectation values of the two Higgs 
fields tan/3 = 15, and the sign of the Higgsino mass term 
fj, > 0. The neutralino lifetime is not defined within the 
model. For this analysis, it is assumed to be sufficiently 
short to yield decays with prompt photons. 

Searches for GMSB SUSY were carried out by collab- 
orations at the CERN LEP collider and at the Fer- 
milab Tevatron collider in both Run I Q and early in 
Run II 0, Q. The initial limits from CDF and DO for 
Run II, based on the SPS 8 model, were combined [9j 
to yield A > 84.6 TeV corresponding to the limit on the 
chargino mass of 209 GeV, at 95% confidence. Comple- 
mentary searches for GMSB SUSY with R-parity viola- 
tion were performed by the HI experiment at HERA [lo| . 



This analysis is an update of that described in Ref. [71] , 
with about a factor of three more data and improved 
photon identification based on: (i) an electromagnetic 
(EM) cluster "pointing" algorithm that predicts the ori- 
gin of a photon with a resolution of about 2 cm along the 
beam axis, thereby eliminating the largest instrumental 
background associated with misreconstruction of the pri- 
mary interaction vertex, and (ii) an improved track veto 
requirement that suppresses sources of background with 
electrons in the final state. We also use an improved like- 
lihood fitter [ll[ to set limits on the scale parameter A. 

The data in this analysis were recorded using sin- 
gle EM triggers with the DO detector [13], the main 
components of which are an inner tracker, liquid- 
argon/uranium calorimeters, and a muon spectrometer. 
The inner tracker consists of silicon microstrip and cen- 
tral scintillating-fiber trackers located in a 2 T supercon- 
ducting solenoidal magnet, providing measurements up 
to pseudorapidities ^ of ]ry] « 3.0 and \ri\ « 1.8, respec- 
tively. The calorimeters are finely segmented and consist 
of a central section (CC) covering \ri\ < 1.2 and two end- 
cap calorimeters extending coverage to \ri\ « 4, all housed 
in separate cryostats [ISJ. The electromagnetic section 
of the calorimeter has four longitudinal layers and trans- 
verse segmentation of 0.1 x 0.1 in rj — (j> space (where (f> 
is the azimuthal angle), except in the third layer, where 
it is 0.05 X 0.05. The central preshower (CPS) system is 
placed between the solenoid and the calorimeter cryostat 
and covers \r]\ < 1.2. The CPS provides precise measure- 
ment of positions of EM showers. The axes of EM show- 
ers are reconstructed by fitting straight lines to shower 
positions measured in the four longitudinal calorimeter 
layers and the CPS (EM "pointing"). The data for this 
study were collected between 2002 and summer 2006, us- 
ing inclusive single EM triggers that are almost 100% 
efficient to select signal data. The integrated luminos- 
ity [li] of the sample is 1100 ± 70 pb-^. 

Photons and electrons are identified based on recon- 
structed EM clusters using calorimetric information and 
further classified into electron and photon candidates, 
based on tracking information. The EM clusters are se- 
lected from calorimeter clusters using the simple cone 
method (of radius 7^ = ^/i^v)^ + (A0)2 = 0.4) by re- 
quiring that (i) at least 90% of the energy is deposited 
in the EM section of the calorimeter, (ii) the calorimeter 
isolation variable I = [EtotiOA) - Eem{0.2)]/ Eem{0.2) 
is less than 0.07, where £'tot(0.4) is the total shower en- 



Pseudorapidity is defined as — log(tan(^)), where 9 is the angle 
between the particle and the proton beam direction. 
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ergy in a cone of radius TZ = 0.4, and -B£;m(0.2) is the 
EM energy in a cone of radius TZ = 0.2, (iii) the trans- 
verse, energy- weighted, width of the EM cluster in the 
third EM calorimeter layer is smaller than 0.04 rad, and 
(iv) the scalar sum of the transverse momenta {j)t) of 
all tracks originating from the primary vertex in an an- 
nulus of 0.05 < TZ < 0.4 around the cluster is less than 
2 GeV. The EM cluster is further defined as an elec- 
tron candidate if it is spatially matched to activity in 
the tracker, and as a photon candidate otherwise. The 
tracker activity can be either a reconstructed track or a 
density of hits in the silicon microstrip and central fiber 
trackers consistent with a track, i.e., an electron. The 
latter requirement allows for increasing electron track- 
matching efficiency, etik, measured in Z — > ee data, from 
(93.0 ± 0.1)% to (98.6 ± 0.1)% by identifying electrons 
with lost tracks due to hard bremsstrahlung and/or in- 
efficiency of the inner trackers. This reduces electron 
backgrounds to photons by a factor of five, while keeping 
the efficiency of anti-track activity requirement high. We 
measure that (91 ± 3)% of photon candidates in Z ^ eej 
data satisfy the anti-track activity requirement. 

Jets are reconstructed using the iterative, midpoint 
cone algorithm [l^ with a cone size of TZ — 0.5. The 
missing transverse energy is determined from the energy 
deposited in the calorimeter for I77I < 4 and is corrected 
for the EM and jet energy scales. 

We select 77 candidates by requiring events to have 
two photon candidates, each with transverse energy 
Et > 25 GeV identified in the CC with \r]\ < 1.1. We 
require that at least one of the photon candidates be 
matched to a CPS cluster, and that the primary vertex 
be consistent with that of the photon candidate (obtained 
from the EM pointing). The accuracy of the determina- 
tion of the photon vertex is measured using photons from 
final state radiation in Z eej data sample and found to 
be 2.3 ± 0.3 cm. The requirement of consistency between 
the photon and primary vertices ensures correct calcu- 
lation of the transverse energies and tracking isolation 
requirements. The accuracy of primary vertex associa- 
tion is studied in GMSB SUSY Monte Carlo simulated 
events, where the primary vertex is identified correctly 
in (98.5 ± 0.1)% of the events while the photon vertex 
matches the primary vertex in (95.8 ± 0.1)%. 

To reduce potential bias in the measurement of 
]^rp from mismeasurement of jet transverse momentum, 
we also require that the jet with the highest Et (if jets 
are present in the event) be separated from the iii 
imuth by no more than 2.5 radians. This selection yields 
2341 events (the 77 sample). 

All instrumental backgrounds arise from standard 
model processes, with either genuine T^^ {Wj, W^-|-jet, 
and tt production) or without inherent T^^p (direct pho- 
ton, multi-jet, and Z ^ ee production). All these back- 
grounds are measured using data. 

The former source always has an electron in the final 



state which is misidentified as a photon. The contribu- 
tion of this background to the l^x distribution in data 
can be estimated using an ej sample (selected by re- 
quiring an electron and a photon candidate and using 
the same kinematical requirements as for the 77 sample) 
scaled by the probability of an electron-photon misiden- 
tification which is measured using Z — *■ ee data. First, 
the l^rp distribution in the 67 sample must be corrected 
for the contribution from events with no real l^x- The 
contribution from Drell-Yan events is taken into account 
by obtaining the l^rp distribution for the ee sample (se- 
lected by requiring two electron candidates and applying 
the same kinematical requirements as for the 77 sample) 
which is dominated by Drell-Yan events. The Drell-Yan 
l^x distribution is further normalized to the number of Z 
boson events in the ej sample (the latter is determined 
by fitting the ej invariant mass spectrum to the Z boson 
mass peak). 

The contribution from the multi-jet processes is esti- 
mated from a data sample (referred to as the QCD sam- 
ple) selected by requiring two EM clusters that (a) sat- 
isfy all the kinematic selection used to select 77 sample 
and (b) satisfy all the photon identification criteria but 
fail the shower shape requirement. The J^^ distribution 
in the QCD sample is normalized to the number of the 
events in the ej sample with T^j^ < 12 GeV after subtrac- 
tion of the Drell-Yan contribution as determined above. 
The expected number of W-f, W+jet, and tt events with 
< 12 GeV is negligible. 

After the Drell-Yan and multi-jet contributions to the 
67 sample are subtracted, the resulting l^rp distribution 
is scaled by (1 — etrk)/etrk, where etrk is the efficiency of 
the track-matching requirement to obtain the estimate of 
distribution for the background with genuine l^x ■ 

The background from events with no inherent T^rp is 
divided into multi-jet events with two real isolated pho- 
tons and events where one or both photons are misiden- 
tified jets. Since the l^x resolution for both sources 
is dominated by the photon energy resolution, the 
distributions for the two sources are very similar. 
However, misidentified jets have a different energy re- 
sponse compared with that of real photons which leads to 
a slight difference in the shapes of the distributions. 
For the real di-photon events, the l^x is assumed to 
have the same shape as that of the Drell-Yan events. 
For misidentified jets, the shape of the l^rp distribution 
is taken from the QCD sample. Relative normaliza- 
tion of the two sources is obtained using a fit to the 
distribution in the 77 sample. We check that the 
fit is not sensitive to possible signal contribution, and 
cross-check with a method that estimates the 77 sam- 
ple purity using the measured shower shape in the CPS. 
The relative fraction of di-photons is (60 ± 20)% and 
this uncertainty is propagated as a systematic uncer- 
tainty for the limit setting. Absolute normalization of 
the T^rp distributions from both sources is determined so 
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n 
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FIG. 1: The l^j, distribution in 77 data with W/Z + -y^ back- 
ground (hatched histogram), instrumental background with 
no genuine J^j,: 77 (sohd black line) and multi-jet (filled 
histogram), and background from processes with genuine 
and a misidentified electron (cross-hatched histogram). 
The expected ]^rp distributions if GMSB SUSY events were 
present are shown as dotted and dashed lines. 



that the number of events with l^j- < 12 GeV matches 
that in the 77 sample. 

The largest physics backgrounds are from Z77 
!/z^77 and Wjj £771/ processes. Contribu- 

tions from these backgrounds are estimated as 0.15 ± 
0.06 and 0.10 ± 0.04 events, respectively, using 
CompHep fl^ Monte Carlo simulation, cross-checked 
with MadGraph fvi\. The contribution of these back- 
grounds to the distribution is taken from Monte Carlo 
simulation, with number of events normalized to the in- 
tegrated luminosity of the data sample. 

The distribution for the 77 sample, with contri- 
butions from physics background {W/Z + 77), and in- 
strumental background with genuine l^rp (processes with 
misidentified electrons) and no inherent l^j- (77 and 
multi-jet) is given in Fig [1] We also illustrate the 
distribution expected from GMSB SUSY for two val- 
ues of A. The number of observed events, as well as 
expected background and signal from GMSB SUSY for 
> 30 GeV and > 60 GeV are given in Table U 

The expected GMSB signal efficiency is estimated from 
Monte Carlo simulation generated for several points on 
the Snowmass Slope (see Table |lTl), covering the neu- 
tralino mass range from 170 GeV to 280 GeV. Al- 
though xtxi &n.d xf X2 processes dominate, we con- 
sider all GMSB SUSY production channels. We used 
IS A JET 7.58 [iil to determine SUSY interaction eigen- 
state masses and couplings. PYTHIA 6.319 [31 is used 
to generate the events after determining the sparticle 
masses, branching fractions and leading order (LO) pro- 
duction cross sections using CTEQ6L1 parton distribu- 
tions [2^. The generated events are processed through 



a full GEANT-based '2l'l detector simulation and the 
same reconstruction code as used for data. The LO signal 
cross sections are scaled to match the next-to-leading or- 
der (NLO) prediction using fc-factor values (see Table lll|. 
extracted from Ref . [22] . 

The systematic error on the expected number of signal 
events comes from the uncertainties in photon identifi- 
cation efficiency (10%), statistics in MC samples (5%), 
track veto requirement (3%), and trigger efficiency (4%). 
These were obtained using Z — » e~^e~ and Z — > e~^e~j 
decays in data and in MC simulation. Variation of par- 
ton distribution functions and uncertainty in the total 
integrated luminosity result in additional 4% and 6.1% 
errors in signal yield respectively. The total uncertainty 
on the background is dominated by statistics. 

As the observed number of events for all values of 
is in good agreement with the standard model pre- 
diction, we conclude that there is no evidence for GMSB 
SUSY in the data. We set limits on the production cross 
section by utilizing a likelihood fitter jH)] that incorpo- 
rates a log-likelihood ratio (LLR) test statistic method. 
This method utilizes binned ]^rp distributions rather than 
a single-bin (fully-integrated) value, and therefore ac- 
counts for the shapes of the distributions, leading to 
greater sensitivity. The value of the confidence level for 
the signal CLg is defined as CLg = CLs+b/CLi,, where 
CLs+b and CLb arc the confidence levels for the sig- 
nal plus background hypothesis and the background-only 
(null) hypothesis, respectively. These confidence levels 
are evaluated by integrating corresponding LLR distri- 
butions populated by simulating outcomes via Poisson 
statistics. Systematic uncertainties are treated as un- 
certainties on the expected numbers of signal and back- 
ground events, not the outcomes of the limit calculations. 
The degrading effects of systematic uncertainties are re- 
duced by introducing a maximum likelihood fit to the 
missing transverse energy distribution. A separate fit is 
performed for both the background-only and signal-plus- 
background hypotheses for each data or pseudo-data dis- 
tribution. 

The limits are shown in Fig. [2] together with expected 
signal cross sections. The observed limits are statistically 
compatible with the expected limits. The observed upper 
limit on the signal cross section is below the prediction 
of the Snowmass Slope model for A < 91.5 TeV, or in 
terms of gaugino masses, m^o < 125 GeV and m-+ < 
229 GeV. These represent the most stringent limits on 
this particular GMSB SUSY model to date. 
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Background events 




Expected 


signal events 


Observed events 




Genuine 


No Physics 


Total 


A = 75 TeV 


A = 90 TeV 




$T > 30 GeV 


0.97±0.12 


9.62±1.12 0.19±0.07 


10.8±1.1 


28.3±4.2 


8.7±1.3 


16 


> 60 GeV 


0.11±0.04 


1.44±0.43 0.08±0.04 


1.6±0.4 


18.1±2.7 


6.4±1.0 


3 



TABLE L Numbers of background events from W^, W + jet, and ti (Genuine ]^t)i no inherent (No -^y), ^77 viyyy 
and VK77 — > ^771/ (Physics) processes; the total number of expected background events; numbers of expected GMSB SUSY 
signal events for two values of A; and the observed numbers of events for > 30 GeV and 60 GeV. Errors are statistical and 
systematic combined. 



, TeV 


m^o, GeV 


171.+ , GeV 


a^°, fb fc- factor 


Efficiency 


70 


93.7 


168.2 


215 


1.21 


0.17 ±0.03 


75 


101.0 


182.3 


148 


1.20 


0.18 ±0.03 


80 


108.5 


198.1 


97.5 


1.19 


0.18 ±0.03 


85 


115.8 


212.0 


65.4 


1.18 


0.19 ±0.03 


90 


123.0 


225.8 


41.8 


1.17 


0.19 ±0.03 


95 


130.2 


239.7 


29.5 


1.16 


0.20 ±0.03 


100 


137.4 


253.4 


20.6 


1.15 


0.20 ±0.03 


105 


144.5 


267.0 


UA 


1.14 


0.18 ±0.03 


110 


151.7 


280.7 


10.3 


1.13 


0.19 ±0.03 



TABLE II: Points on the GMSB Snowmass Slope model: 
neutralino and chargino masses, cross sections predicted by 
PYTHIA, fc-factors, and reconstruction efficiencies with total 
uncertainty. 
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FIG. 2: Predicted cross section for the Snowmass Slope model 
versus A. The observed and expected 95% C.L. limits are 
shown in solid and dash-dotted lines, respectively. 
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